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ABSTRACT 

i We have derived star formation rates (SFRs), gas-phase oxygen abundances and ef- 

' fective dust absorption optical depths for a sample of galaxies drawn from the Stromlo- 

, APM redshift survey using the new Chariot & Longhetti (2001; CLOl) models, which 

provide a physically consistent description of the effects of stars, gas and dust on the 
integrated spectra of galaxies. Our sample consists of 705 galaxies with measurements 
of the fluxes and equivalent widths of Ha, [On], and one or both of [Nil] and [Sii]. 
For a subset of the galaxies, 60 and 100 fim IRAS fluxes are available. We compare 
fH I the star formation rates derived using the models with those derived using standard 

Q^i estimators based on the Ha, the [O ii] and the far- infrared luminosities of the galaxies. 

I ' The CLOl SFR estimates agree well with those derived from the IRAS fluxes, but are 

O I typically a factor of ~ 3 higher than those derived from the Ha or the [On] fluxes, 

even after the usual mean attenuation correction of = 1 mag is applied to the 
^ ' data. We show that the reason for this discrepancy is that the standard Ha estimator 

neglects the absorption of ionizing photons by dust in H II regions and the contamina- 
tion of Ha emission by stellar absorption. We also use our sample to study variations 
in star formation and metallicity as a function of galaxy absolute bj magnitude. For 
this sample, the star formation rate per unit bj luminosity is independent of magni- 
' tude. The gas-phase oxygen abundance does increase with bj luminosity, although the 

scatter in metallicity at fixed magnitude is large. 

Key words: galaxies: general - galaxies: ISM - galaxies: stellar content. 



1 INTRODUCTION 

Most large spectroscopic galaxy surveys have been carried 
out with the aim of determining the luminosity function 
or of studying the three-dimensional clustering properties 
of galaxies. The spectra are used to measure redshifts and 
so to infer the distances to galaxies. However, the spectra 
also contain a wealth of additional information about the 
physical properties of galaxies. The continuum and the ab- 
sorption lines provide information about the stellar content. 



* Email: charlot@iap.fr 
© 1994 RAS 



while the nebular emission lines provide a measure of the 
star formation rate (SFR) and the interstellar metallicity. 

Kennicutt (1992a) published an atlas of high-resolution, 
high signal-to-noise, flux-calibrated spectra of a small sam- 
ple of nearby galaxies of different Hubble types. The spec- 
tra obtained in redshift surveys are typically of much lower 
quality than the spectra in the Kennicutt atlas, but they 
nevertheless provide a means of studying the star formation 
rates and metallicities of a magnitude-limited sample of nor- 
mal galaxies and of characterizing trends in these properties 
as a function of galaxy luminosity and environment. As a 
result, there have been a large number of papers analyzing 
the spectral properties of redshift survey galaxies. 

Many studies focus on SFR estimates based on a single 
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strong emission line, such as [O ll]A3727 (e.g. Hashimoto et 
al 1998; Christlein 2000; Baldi, Bardelh & Zucca 2001). One 
of the advantages of [O ll] is that it is easily observed in high- 
redshift galaxies and can thus be used to study evolutionary 
trends in star formation (e.g. Ellis et al 1996; Cowie et al 
1996; Hammer et al 1997; Balogh et al 1998). Other studies 
focus on additional lines, in particular the Ha line, which is 
believed to be a more robust tracer of star formation than 
[O ll] (e.g. Gallego et al 1995; Tresse & Maddox 1998; Tresse 
et al 1999; Barton, Geller & Kenyon 2000). 

Kennicutt (1998; hereafter K98) reviewed a number of 
traditional diagnostic methods for inferring star formation 
rates from galaxy spectra. The two methods that are appli- 
cable to redshift survey data are based on the Ha recombina- 
tion line and the [O ll] forbidden-line doublet. As discussed 
by K98, the derived SFRs are subject to errors because of 
attenuation by dust. In the case of [O ll], variations in excita- 
tion, metallicity and diffuse gas fraction will also contribute 
to the uncertainty in the derived SFRs. Both methods ignore 
the absorption of ionizing photons by dust in Hll regions, 
which could play an important role (Petrosian, Silk & Field 
1972; Mathis 1986; Chariot & FaU 2000). 

Recently, Chariot & Longhetti (2001, hereafter CLOl) 
quantified the uncertainties in these SFR estimators. They 
constructed physically realistic models describing the propa- 
gation of the photons emitted by stars through the interstel- 
lar medium (ISM) and their absorption by dust in galaxies. 
These models combine the latest version of the Bruzual & 
Chariot (1993) population synthesis code with the Ferland 
(1996) photoionization code for modelling nebular emission. 
Stars are assumed to form in interstellar birth clouds (i.e. 
giant molecular clouds). After 10^ yr, young stars disrupt 
their birth clouds and migrate into the ambient ISM. The 
line emission from a galaxy arises both from the H ll regions 
surrounding young stars and from the diffuse gas ionized by 
photons that have leaked into the ambient ISM. In the CLOl 
model, the depletion of heavy elements onto dust grains, the 
absorption of ionizing photons by dust, and the contamina- 
tion of Balmer emission lines by stellar absorption are all 
included in a self-consistent way. The absorption of photons 
from H II regions and from older stars is described by a dust 
model developed by Chariot & Fall (2000). In this prescrip- 
tion, line and continuum photons are attenuated differently 
because of the finite lifetimes of the clouds in which stars 
form. 

CLOl showed that in models that include absorption 
by dust and that match the observed ionized-gas properties 
of local galaxies, SFR estimates based purely on the Ha or 
[O ll] luminosities of galaxies can be in error by more than 
an order of magnitude. On the other hand, with the help of 
other lines such as H/3, [Olll], [Nil] and [Sll], these errors 
can be substantially reduced. CLOl presented a number of 
new estimators of the SFRs, gas-phase oxygen abundances 
and dust absorption optical depths of galaxies based on com- 
binations of these lines. 

This paper presents the first application of the CLOl 
model to a magnitude-limited sample of galaxies drawn from 
the Stromlo-APM (SAPM) redshift survey (Loveday et al 
1996). Tresse et al (1999) have measured the fluxes and 
equivalent widths of Ha, [On], [Nil] and [Sii] for the galax- 
ies in this survey. We compare the SFRs derived using stan- 
dard estimators based on the Ha and [O ll] lines to those 



derived using the CLOl models, which make use of all four 
lines. We show that the CLOl SFR estimates are typically 
a factor of ~ 3 higher than those derived using just Ha or 
[On], even after the usual mean attenuation correction of 
Aua ~ 1 mag is applied to the data. We show that the 
reason for this discrepancy is that the standard Ha estima- 
tor neglects the absorption of ionizing photons by dust in 
H II regions and the contamination of Ha emission by stellar 
absorption. 

Observations from the Infrared Astronomical Satellite 
(IRAS) are available for 149 galaxies in the sample and pro- 
vide an independent check on the derived star formation 
rates. We find that the CLOl SFR estimates agree surpris- 
ingly well with those based on the 60 and 100 /im IRAS 
fluxes. 

Within the subset of the Stromlo-APM survey which we 
study, we find that the star formation rate per unit bj lu- 
minosity is independent of galaxy absolute magnitude. We 
derive gas-phase oxygen abundances and dust absorption 
optical depths for this sample. Bright, L, galaxies have 0/H 
values close to solar, while fainter galaxies are found to be 
less metal-rich on average. The distribution of Ha attenu- 
ations is similar to that found for other samples of nearby 
star- forming galaxies. 

We describe the properties of the SAPM galaxy sample 
in Section 2 below. In Section 3, we illustrate the incon- 
sistencies obtained when estimating star formation rates of 
galaxies using different traditional diagnostic methods (Ha, 
[On], UV and FIR luminosities). Readers familiar with the 
uncertainties affecting star formation measures in galaxies 
may skip at least the first part of this section. In Section 4, 
we describe how the new models developed by CLOl can help 
us better constrain the star formation, metallicity and dust 
properties of SAPM galaxies using a combination of differ- 
ent spectral diagnostics. In Section 5, we present the results 
of this analysis and identify the important weaknesses af- 
fecting the traditional Ha estimator. Our conclusions are 
summarized in Section 6. 



2 DESCRIPTION OF THE SAMPLE 

Our sample is drawn from the the Stromlo-APM (SAPM) 
redshift survey (Loveday et al. 1996). The SAPM survey 
consists of 1797 galaxies brighter than bj = 17.15 cover- 
ing 4300 deg^ of the south galactic cap. The mean redshift 
of the galaxies in the survey is (z) = 0.051. The spectra 
were taken with the Dual-Beam Spectrograph of the Aus- 
tralian National University 2.3-m telescope at Siding Spring. 
The wavelength range is 3700-5000 A in the blue and 6300- 
7600 A in the red, and the dispersion is ~ 1 A/pixel. The 
fibre-coupling of each beam to two CCD chips introduced ad- 
ditional small gaps in the wavelength coverage in the regions 
4360-4370 A and 7000-7020 A. The spectra were taken with 
8 arcsec slits and the spectral resolution was 5 A (FWHM). 

Tresse et al. (1999) measured the 4000 A breaks and 
emission-line properties of 1671 galaxies with 15 < foj < 
17.15 in the sample (the other 126 galaxies in the original 
sample were either brighter than bj = 15, or their spectra 
were blue-shifted, or they had published redshifts). They 
measured the fiuxes of all lines with emission equivalent 
widths larger than 2 A, corresponding to a 3a detection 
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limit. Because of the wavelength gaps in the spectra, differ- 
ent lines were measured for galaxies at different redshifts. 

We first exclude 79 galaxies with [N n]A6583/HQ: > 
0.63, which are likely to be AGN. The fluxes of Ha, 
[O n]A3727, [N ii] A6583 and [S ii] AA6717, 6731 could be mea- 
sured for 416 galaxies. The fluxes of subsets of these lines are 
available for 588 additional galaxies. Of these, 289 galaxies 
have measurements of Ha, [On], and either [Nil] or [Sll]. 
Finally, 588 galaxies have no fluxes in emission with equiv- 
alent widths larger than 2 A. We note that since H/3 and 
[Oni]A5007 fall in between the blue and red beams of the 
spectrograph for most objects, they are generally not ob- 
served. In this paper, we focus on the subsample of 705 non- 
Seyfert galaxies for which Hq, [O ii], and at least one of [N ll] 
or [S ii] were measured. 

The absolute flux calibration of the original SAPM 
spectra is only accurate to ~ 20 per cent because of the 
limited number of standard stars observed each night (Sin- 
gleton 2001). However, the relative flux calibration, which 
is more critical for our analysis of emission line ratios and 
equivalent widths is accurate to ~ 5 per cent (Tresse et al. 
1999). The 8 arcsec-wide and 7 arcmin-long slit used in the 
observations samples typically 45 per cent of the rest-frame 
projected area of a galaxy at a b,j surface brightness level of 
25 mag arcsec"^. To correct for this effect. Singleton (2001) 
computed the flux in the bj band in each spectrum and de- 
rived an aperture correction for each galaxy using the ratio 
of the spectral flux to the total bj band flux of the galaxy. 

The above aperture correction should be accurate so 
long as the line emission and the blue light are correlated 
over the entire galaxy. One worry is that redshift surveys 
that use limited spectroscopic apertures may be biased as a 
function of morphological type, because spectroscopic prop- 
erties vary significantly from the inner (bulge-dominated) to 
the outer (disc-dominated) regions of galaxies. Kochanek, 
Pahre & Falco (2001) have shown, however, that the long- 
slit SAPM spectra sample the overall light distribution of 
a galaxy far better than the small-aperture (e.g. 2-3 arcsec 
fibre) spectra of other nearby spectroscopic surveys. They 
conclude that the spectral classification of SAPM galax- 
ies should not suffer from aperture bias. Also, Tresse et 
al. (1999) showed that the average properties of galaxies 
in the SAPM sample are similar, in terms of line ratios 
(e.g. [Nii]/Hq: and [Oll]/Ha), to those in the integrated- 
spectrum sample of Kennicutt (1992b). As found by Ken- 
nicutt (1992b), the main effect of disc undersampling is to 
reduce the strengths of the various emission lines in roughly 
equal proportion. Thus, the aperture corrections derived by 
Singleton (2001) should be appropriate. We expect, there- 
fore, that the spectral analysis of the SAPM galaxy sample 
presented below does not suffer from any significant aperture 
bias. 

Observations from IRAS are also available for 279 galax- 
ies in the sample (Singleton 2001). Only 149 of these belong 
to our subsample of 705 non-Seyfert galaxies with measure- 
ments of Hq, [Oil], and one or both of [Nil] and [Sll]. 

In summary, we consider two samples in this paper: (1) 
a restricted sample of 149 non-Seyfert galaxies with IRAS 
detections, for which fluxes and equivalent widths of Ha, 
[O ii], and one or both of [N ii] and [S ii] are available. We use 
this sample to compare star formation rates derived using 
different estimators based on the Hq, [O ii] and far-infrared 



(FIR) luminosities of galaxies; (2) a larger sample of 705 
non-Seyfert galaxies for which ffuxes and equivalent widths 
of Ha, [O ii], and one or both of [N ll] and [S ll] are available. 
We use this sample to study trends in star formation rate 
and metallicity as a function of galaxy absolute bj luminos- 
ity. 

In this paper, we adopt a Hubble constant Ho = 
70 kms"^Mpc"^ 



3 STAR FORMATION RATES DERIVED 
USING STANDARD ESTIMATORS 

In this section, we adopt the formulae given in K98 to con- 
vert Ha, [O ii], UV and FIR luminosities into star formation 
rates, and we check the results for consistency. 

3.1 Standard SFR estimators 

The K98 calibration of the Ha luminosity emitted per unit 
SFR is derived by using standard population synthesis mod- 
els to calculate the ionizing radiation produced by young 
stars and by applying dust-free case B recombination to this. 
For solar abundances and a Salpeter IMF (0.1-100 M©), K98 
derives the following transformation between Ha luminosity 
and SFR: 

SFRhc (Moyr-') = 7.9 x IO-^'Lhc (ergs"'). (1) 

This calibration does not include the effects of at- 
tenuation by dust. Kennicutt (1983) and Niklas, Klein & 
Wielebinski (1997) have used observations of the integrated 
Ha and thermal radio fluxes of nearby spiral galaxies to 
derive a mean attenuation Ana ~ 0.8 — 1.1 mag for these 
objects. Unless otherwise specifled, we apply a 'standard' 
attenuation correction Aua = 1 mag to all galaxies in the 
sample. We emphasize that this corrects for only part of the 
effects of dust on the Ha line, namely the absorption by 
dust of Ha photons outside the Hll regions in which they 
are produced. The correction neglects the diminution of the 
Ha line caused by the absorption of ionizing photons by 
dust inside Hll regions (Petrosian, Silk & Field 1972; Mathis 
1986; Chariot & Fall 2000). More reflned attenuation correc- 
tions based on the departure of the observed Ha/H/3 ratio of 
each galaxy from dust-free case B recombination would suf- 
fer from the same limitation (we recall that H/3 luminosities 
are not available for SAPM galaxies). 

The luminosities of forbidden lines such as [O ll] are not 
directly coupled to the ionization rate, and their excitation 
is known to be sensitive to the abundance and the ionization 
state of the gas. However, it is claimed that the excitation of 
[O ii] is sufficiently well behaved in observed galaxies that it 
can be calibrated empirically through Ha as a quantitative 
SFR tracer (Gallagher et al 1989; Kennicutt 1992b). K98 
gives the formula 

SFR[o„] (Moyr-i) = (1.4 ± 0.4) x 10-"L[o„] (ergs"^). (2) 

It is worth noting that, because this estimator is based on an 
empirical calibration, the relative attenuation of the Ha and 
[O ii] lines is automatically accounted for. We also apply a 
standard attenuation correction Ana ~ 1 mag to the [O ii]- 
derived SFR estimates. 
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The non-ionizing ultraviolet (UV) radiation from young 
stars is another standard SFR estimator, especially in high- 
redshift galaxies where optical emission lines are not always 
accessible. K98 gives the following transformation between 
UV luminosity L^,uv and SFR: 

SFRuv(M0yr-') = 1.4 x 10"^**L,,uv (erg s-^Hz"^). (3) 

This formula includes no dust correction and assumes that 
the SFR has remained constant over timescales that are long 
compared to the lifetimes of the dominant UV emitting pop- 
ulation (< 10* yr). For the choice of a Salpeter (0.1-100 Mq) 
IMF, the UV luminosity is practically independent of wave- 
length. In Section 3.2 below, we address the uncertainty in 
the above relation that arises as a result of the absorption 
of UV radiation by dust. 

The FIR emission should provide an excellent mea- 
sure of the SFR in dusty, circumnuclear starbursts. As dis- 
cussed by K98, in the optically thick limit, it is sufficient to 
model the bolometric luminosity of the stellar population, 
the greatest uncertainty being the assumed age of the star- 
burst. K98 quote the following relation for continuous bursts 
of age 10-100 Myr: 

SFRfir (M0 yr-') = 4.5 x IO-^^Lfir (ergs"') , (4) 

where Lfir refers to the luminosity integrated over the full 
mid-infrared to submillimeter spectrum (8-1000 /im). 

In more normal star-forming galaxies, the relation be- 
tween SFR and FIR emission is more complex. As discussed 
by K98, the contribution to dust heating from older stars 
will lower the effective coefficient in Equation (^, whereas 
the lower optical depth of the dust will tend to increase 
the coefficient. Buat & Xu (1996) computed the total FIR 
luminosity emitted per unit SFR in spiral galaxies of type 
Sb to Irr using ultraviolet and IRAS observations of 152 
nearby disc galaxies. For this sample, they derived a mean 
coefficient a factor ~ 2 larger than that given in Equation 
. Furthermore, Misiriotis et al. (2001) have carried out 
a detailed study of the optical to FIR spectral energy dis- 
tributions of nearby, edge-on galaxies of type Sb-Sbc. Their 
data also indicate that naive application of Equation (^) to 
these galaxies will underestimate the SFRs by a factor of 
~ 2. [Note that despite these considerations, we will always 
adopt the expression in Equation above when computing 
SFRfir.] 

3.2 Comparison of SFRs derived using different 
standard estimators 

The SFR estimators above are valid in opposite limits for the 
transfer of line and continuum radiation through the ISM of 
galaxies. Equation (^ assumes that the galaxies are trans- 
parent. Standard attenuation corrections are at best very 
approximate remedies. Equation ^ simply rewrites Equa- 
tion (|^) using an empirically derived proportionality con- 
stant, and so applies in the same limit (although with sig- 
nificant additional uncertainties). Equation (^), which also 
applies to transparent galaxies, traces the emission from 
stars that live longer than those dominating the ionizing 
radiation. Equation (^ , on the other hand, is valid only for 
opaque galaxies. It would thus be surprising if the four esti- 
mators agreed, even for galaxies with perfect measurements 
of Lua, L[ou], ii^.uv and Lfir. 



To compare the star formation rates obtained using 
these different estimators, we first estimate the total FIR 
luminosity Lfir for the 149 galaxies in our restricted SAPM 
sample with IRAS flux densities at 60 and 100 fim. We 
adopt the following standard procedure (see for example 
Meurer et al. 1999). We first calculate the quantity Ffir = 
1.26 X 10"" [2.58/^ (60 /xm) + /^(lOO/im)] ergcm'^s"' de- 
fined by Helou et al. (1988) and then compute the total far- 
infrared flux from Fpin and (60 /im) / /i^ (100 fim) using the 
relation appropriate for dust with a single temperature and 
an emissivity proportional to frequency f. This prescription 
has been shown work well for starburst galaxies (see Fig. 2 
of Meurer et al. 1999) . It is expected to a good approxima- 
tion for all galaxies in which the infrared emission peaks at 
wavelengths between 60 /im and 100 /im (Helou et al. 1988). 
The uncertainties that arise from the conversion of IRAS 
fluxes to total FIR flux are likely to be modest in compari- 
son to the uncertainties in the conversion from Lfir to star 
formation rate. 

Fig. 1 shows the logarithmic difference between SFRhq, 
SFR[oii] and SFRfir for the 149 SAPM galaxies in this 
sample as a function of absolute rest-frame bj magnitude. 
In each panel, we indicate the median A log (SFR) offset for 
galaxies with M{bj) < —19 (representing 80 per cent of 
the sample), along with the associated mean offset and the 
rms scatter (in parentheses). For both Hq and [On], we 
have included the 'recommended' mean attenuation correc- 
tion Ana ~ 1 mag. As expected, the three estimators yield 
inconsistent results. 

Fig. la is consistent with the recent result by Jansen, 
Franx & Fabricant (2001; see also Carter et al. 2001) that 
bright galaxies tend to have lower [On] /Ha ratios than less 
luminous galaxies. We find that, for M{bj) < —19, [Oll]- 
derived SFR estimates are typically 40 per cent smaller than 
Ha-derived ones, while at fainter magnitudes, the two es- 
timators give similar results. At fixed magnitude, the rms 
scatter in SFR[oii]/SFRHa is a factor of ~ 1.7, although 
discrepancies between Ha- and [O n]-derived SFR estimates 
can reach a factor of 7. The scatter in this plot and the weak 
trend with magnitude can probably be explained by varia- 
tions in the effective gas parameters (ionization, metallicity, 
dust content) of the galaxies. We note that M{bj) ~ —19.5 
corresponds to the mean magnitude of the galaxies which 
K98 used to calibrate [O ll] as a SFR estimator. 

The differences between the SFRs estimated using op- 
tical lines and those estimated using the total FIR luminos- 
ity are much more severe. Figs. Ib-lc show that Ha- and 
[O n]-derived SFR estimates are typically a factor of ~ 3 
lower than FIR-derived ones, even though we have applied 
an attenuation of 1 mag at Ha to the data. We note that, 
as discussed above, the K98 formula given in Equation (^ 
may well underestimate the true star formation rate in nor- 
mal star-forming galaxies. As we shall see in Section 5, the 
discrepancy in Figs. Ib-lc is likely to arise from the conver- 
sion of optical emission-line luminosities into star formation 
rates. 

Ultraviolet luminosities are not available for the SAPM 
galaxies. To compare SFRh^ and SFRfir to the star for- 
mation rates derived using UV luminosities, we have an- 
alyzed 37 UV-selected, normal star-forming galaxies from 
the sample of Bell & Kennicutt (2001) and 18 UV-selected 
starburst galaxies from the sample of Meurer et al. (1999). 
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We selected in the original samples all the galaxies with Ha, 
UV and IRAS 60 /xm and 100 /xm fluxes. The published rest- 
frame UV fluxes correspond to the mean flux density over a 
passband with a central wavelength of Auv ~ 1600 A and a 
width of a few hundred angstrom. 

Fig. 2 shows the logarithmic difference between SFRuv, 
SFRhq and SFRpir for the galaxies in these samples as a 
function of absolute rest-frame B magnitude. For consis- 
tency, the Ha-derived SFRs in Fig. 2a are not corrected for 
the attenuation of Ha photons by dust. Starburst galaxies 
follow roughly the same trends as normal star-forming galax- 
ies in these diagrams, although they have somewhat lower 
UV-derived SFRs at fixed B-band magnitude. This may be 
explained if UV photons are more strongly attenuated by 
dust in starburst galaxies, or if long-lived stars with ages 
10^-10* yr contribute a larger fraction of the UV-luminosity 
in normal star-forming galaxies. Fig. 2b is also consistent 
with an increase of attenuation by dust with increasing lu- 
minosity in nearby starburst galaxies (e.g. Heckman et al. 
1998). Overall, the similarity between the properties of star- 
burst and normal star-forming galaxies in Fig. 2 suggests 
that the starbursts do not strongly influence the global ISM 
parameters of galaxies. 

The median difference between the UV- and Ha-derived 
SFRs, both uncorrected for dust, is small, i.e. less than a fac- 
tor of 2. This was already noted by Bell & Kennicutt (2001) 
for their sample. This means that even though the overall 
attenuation by dust can be large (as Fig. 2b demonstrates), 
it is similar for Ha and ultraviolet photons. Ultraviolet pho- 
tons are expected to be more attenuated than Ha photons 
for a normal extinction curve. On the other hand, the stars 
dominating the Ha emission are more obscured than those 
dominating the UV emission, because they have lifetimes 
(;^ 3 X 10^ yr) shorter than the typical timescale for the dis- 
persal of clouds in which they form [the dispersal timescale 
of giant molecular clouds in the Milky Way is estimated 
to be about lO'^yr; Blitz & Shu (1980)]. Similar results are 
found for nearby starburst galaxies, where the attenuation 
inferred from the Ha/H/3 ratio is typically higher thmi that 
inferred from the ultraviolet and optical spectral continuum 
(e.g., Fanelh, O'Connell & Thuan 1988; Calzetti, Kinney & 
Storchi-Bergmann 1994; Calzetti 1997; Poggianti et al. 1999; 
see Chariot & Fall 2000) . In summary, our results show that 
even when the attenuation is strong, the ratio of the emer- 
gent ultraviolet to Ha luminosities can be close to that in 
the dust-free case, leading to small apparent discrepancies 
between the uncorrected UV- and Ha-derived SFRs. 



4 THE CHARLOT & LONGHETTI MODEL 

In the previous section, we investigated a number of stan- 
dard SFR estimators based on the Ha, [O ll], UV and FIR lu- 
minosities of galaxies. The different estimators give inconsis- 
tent results. These results highlight the need for a physically 
realistic model to interpret the combined signatures of stars, 
gas, and dust in the integrated spectra of galaxies. Such a 
model was recently developed by Chariot & Longhetti (2001; 
hereafter CLOl). 



4.1 Description of the model 



The CLOl model is based on a combination of the latest ver- 
sion of the Bruzual & Chariot (1993) population synthesis 
code and the Ferland (1996, version C90.04) photoionization 
code. Stars are assumed to form in interstellar 'birth clouds' 
(i.e. giant molecular clouds). After 10^ yr, young stars are 
assumed to disrupt their birth clouds and migrate into the 
'ambient ISM'. The line emission from a galaxy arises both 
from the H ll regions surrounding young stars and from the 
diffuse gas ionized by photons that have leaked into the 
ambient ISM. In the CLOl model, these different gas com- 
ponents are combined and are described in terms of 'effec- 
tive' gas parameters. This is motivated by the fact that the 
optical-line ratios in the integrated spectra of nearby spiral 
galaxies are similar to those in the spectra of individual H ll 
regions (Kobulnicky et al. 1999; see CLOl for more detail). 
The spectral evolution of a galaxy is obtained by convolv- 
ing the spectral evolution of a single stellar generation with 
the SFR function. Only stars younger than about 3 x 10® yr 
contribute significantly to the line emission. 



The main adjustable parameters that affect the emis- 
sion from the stars in a galaxy are the IMF, the SFR and 
the stellar metallicity Z. The main parameters affecting the 
emission from the photoionized gas are the effective inter- 
stellar metallicity Z, the zero-age effective ionization param- 
eter {Uo), the effective dust-to- heavy element ratio and 
the effective gas density nn- Z, which was originally called 
'effective gas metallicity' in CLOl, includes all the heavy ele- 
ments in the gas and solid (i.e. dust) phases. It is taken to be 
the same as the metallicity of the ionizing stars, i.e. Z = Z. 
(Uo) describes the typical ionization state of photoionized 
gas in the galaxy. is defined as the mass fraction of heavy 
elements locked into dust grains in the ionized gas. The other 
main parameter of the model is the age of the galaxy t. In 
all models, we adopt a Salpeter IMF truncated at 0.1 and 
100 Mq and we fix hn ~ 30 cm"'', a value that is typical 
of the gas density in Galactic and extragalactic H ll regions. 
Changes in fin have a minor effect on the conversion between 
line luminosities and SFR (see CLOl). CLOl calibrated the 
emission-line properties of their model using the observed 
[Olll]/H/3, [Oll]/[Olll], [Sll]/Ha and [Nll]/[Sll] ratios of a 
representative sample of 92 nearby spiral and irregular, star- 
burst and Hll galaxies. 



The CLOl model also includes the prescription for the 
absorption of photons from H ll regions and from older stars 
that was introduced by Chariot & Fall (2000). In this pre- 
scription, line and continuum photons are attenuated dif- 
ferently because of the finite lifetimes of the stellar birth 
clouds. The effective absorption curve is given by the fol- 
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lowing formula|j 

J fv (A/5500 A) ' , fort'<10^yr, 
' l^fv (A/5500 A) fort'>10Vr, ^ ' 

where t' is the age of any single stellar generation. The wave- 
length dependence of the effective absorption curve t\ is 
constrained by the observed relation between the ratio of 
FIR to UV luminosities and the UV spectral slope for star- 
burst galaxies. The age 10^ yr corresponds to the typical 
lifetime of a stellar birth cloud (i.e. giant molecular cloud). 
The adjustable parameter /x defines the fraction of the to- 
tal dust absorption optical depth of the galaxy contributed 
by the ambient ISM. Chariot & Fall (2000) showed that 
pi ~ 1/3 is needed to reproduce the observed mean relation 
between Ha/H/9 ratio and ultraviolet spectral slope (defined 
near 1600 A) for starburst galaxies. The scatter in the re- 
lation is reproduced if ^ ranges from values of ~ 1 (i.e., all 
the radiation from the Hll regions leaks into the ambient 
ISM) to ~ 1/5 in different galaxies. In this prescription, the 
attenuation of the Ha line radiation, which is produced by 
stars younger than ~ 3 x 10^ yr, can be expressed simply 
as Akc = 1.086 TV (AHa/5500A)-° '' (Ahq = 6563 A). The 
attenuation Av of the 5500 A continuum radiation, which 
is produced by stars both younger and older than 10^ yr, 
depends on the star formation history. 

The different parameters of the CLOl model each have 
a specific influence on the various line and continuum spec- 
tral features of galaxies. This allows one to derive 'physi- 
cal' quantities such as star formation rates, gas-phase oxy- 
gen abundances and effective dust absorption optical depths 
from the optical spectra. The constraints will be more accu- 
rate if more spectral features are available. 



4.2 Application to the SAPM survey 

For the SAPM galaxies, the spectral features available to 
constrain physical quantities are the fluxes and equivalent 
widths of Ha, [On], and one or both of [Nil] and [Sll]. 
For a subset of the galaxies, the total FIR luminosity is 
also available. The parameters we wish to constrain are the 
following. For the gas, they are the interstellar metallicity, 
the ionization parameter, the dust-to-heavy element ratio, 
the effective dust absorption optical depth of the 'neutral' 
ISM (i.e. the birth clouds plus the ambient ISM), and the 
fraction of this contributed by the ambient ISM. For the 
stars, they are the age of the galaxy and the star formation 
history (see Table 1). 

For the SAPM sample, the primary constraint on the 
interstellar metallicity comes from the [Nll]/[Sll] ratio, if 
both metal lines are available. If only one line is available, 

t In Chariot & Fall (2000), Ty is used to denote the total ef- 
fective y-band optical depth, i.e. the sum of the effective optical 
depths in the H n regions, the surrounding H I regions of the birth 
clouds and the ambient ISM (fy = r^" + t™ + f|?^). In the 
CLOl model, the absorption in the Hll regions is computed by 
the photoionization code, and the Chariot & Fall (2000) model 
is used with t™^ = to compute the subsequent absorption in 
the surrounding H I regions of the birth clouds and the ambient 
ISM. Thus, in CLOl and here, Ty represents the effective V-band 
optical depth in the 'neutral' ISM, i.e. ry 



-HI . 



-ISM 

w ■ 



the metallicity will be set by [Nll]/Ha or by [Sll]/Ha. The 
[Sll]/Ha ratio also constrains the dust-to-heavy element ra- 
tio ^d, and the [On]/[Sll] ratio constraints the ionization 
parameter {Uo). The effective absorption optical depth of 
the dust Ty is constrained primarily by the [On] /Ha ra- 
tio. As can be seen, there are at most 3 independent line 
ratios available to constrain 4 gas parameters, so the result- 
ing errors on some of these parameters will be large. The 
constraints on the effective dust absorption optical depth, 
for example, would be far tighter if the Ha/H/3 ratio were 
available. In fact, for the subset of galaxies with far-infrared 
luminosities, the Lhc/Lfib. ratio greatly improves the con- 
straints on Ty. The parameters jj, and t and the history of 
star formation SFR(t) are constrained only weakly by the 
measured equivalent widths of the various lines, the bj ab- 
solute magnitude and the 4000 A break. This is why we do 
not show any results for the parameters pertaining to long- 
lived stars in this paper. The specific influence of the various 
parameters on observables allows one to construct 'optimal' 
estimators of the star formation rate, the gas-phase oxygen 
abundance and the dust absorption optical depth based on 
observed spectral features. This was the approach adopted 
by CLOl (their table 2).^ More accurate constraints and bet- 
ter estimates of the errors in the derived quantities may be 
obtained if the models are used to fit the spectral properties 
of each individual galaxy. This is the procedure we adopt in 
this paper. 

Following CLOl, we first select a set of models that 
span the full range of observed optical-line ratios, ultraviolet 
and far-infrared luminosities of nearby star-forming galax- 
ies. The parameters of these models are listed in Table |l| 
They are the same as those listed in table 1 of CLOl, except 
that we have now extended the ranges of metallicities and of 
effective dust absorption optical depths. We construct a grid 
of 1.6 X 10* models by systematically varying the parameters 
listed in Table |l|. In practice, this is done by interpolating 
from the coarse grid of 7.1 x 10^ models computed by CLOl 
(with extensions in Z and Ty). For each galaxy in our sam- 
ple, we evaluate the goodness of fit of all the models in 
the grid. The SFR, O/H and Ty of the mininum-x^ model 
define our best estimates of the star formation rate, the gas- 
phase oxygen abundance and the effective optical depth of 
the dust in the galaxy. The maximum ranges in SFR, O/H 
and Ty among the models within Ax^ = 1 of the minimum 
provide estimates of the uncertainties in these quantities. 



5 RESULTS 

We now present the star formation rates, gas-phase oxygen 
abundances and dust absorption optical depths we derive 
for the SAPM sample using the CLOl model described in 
Section 4. 

In Fig. 3, we present results for the 149 SAPM galaxies 
with fluxes and equivalent widths of Ha, [On], and one or 
both of [N ll] and [S ll], as well as IRAS 60 and 100 ^m fluxes. 



t We note that the CLOl model reproduces equation (^ to within 
a few percent if = Ty = (i.e. in the case of a completely 
transparent galaxy), and if the effects of stellar Ho absorption 
are ignored. Equation (m is recovered if fiTy 3> 1. 
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Table 1. Parameters of the models reproducing the observed in- 
tegrated spectral properties of nearby star-forming galaxies. Two 
timescales are adopted for the exponentially declining star forma- 
tion rate SFR{t): 0.1 and 6.0 Gyr. All models have nn = 30 cm"^ 
and a Salpeter IMF truncated at 0.1 and 100 Mq. 



Parameter 


Range 


(Uo) 




Z/Zq 


0.2-4.0 


id 


0.1-0.5 


TV 


0.01-4.0 




0.2-1.0 


SFR(t) 


constant, ( 


t/yr 


lOMoio 



In this plot, we do not use the IRAS fluxes to constrain the 
SFRs. We compare the SFRs derived using all available op- 
tical lines with those derived from the Ha, [O ll] or FIR lu- 
minosities alone using the K98 formulae given in Equations 
(0), (P) and (^). The filled squares represent galaxies with 
measurements of all four optical lines, while the open circles 
and the open squares are galaxies for which Ifa, [On], and 
either [N ll] or [S ll] were measured. The error bars indicate 
the uncertainties in the SFRs derived using the CLOl mod- 
els, as discussed in Section 4.2. Once again we have applied 
a standard attenuation Aua = 1 mag to the SFRs estimated 
using the K98 formulae for Ha and [O ll]. 

We find that the SFRs derived using Ha are lower 
than those obtained using the CLOl model. The median dis- 
crepancy tends to increase with bj absolute luminosity and 
amounts to a factor of ~ 3 at M{bj) < — 19 (Fig. 3a). The 
discrepancy is similar, but the scatter is larger for the SFRs 
derived using [On] (Fig. 3b). Interestingly, the CLOl SFRs 
are in much better agreement with those derived from the 
FIR luminosities, even though we dtd not use these measure- 
ments to constrain the models (Fig. 3c). 

Fig. 4 shows what happens if we add the requirement 
that models reproduce the observed FIR luminosities of the 
galaxies. The uncertainties are much smaller in this case 
because fewer models fit the data. The median error on the 
SFR decreases from a factor of ~ 1.8 in Fig. 3 to a factor 
of ~ 1.3 in Fig. 4. The SFRs derived using the CLOl model 
now agree surprisingly well with those derived using the K98 
FIR-luminosity estimator, with a scatter of only 50 per cent. 

Why are the CLOl SFR estimates so different from the 
SFRs derived using the K98 Hq formula plus a standard 
attenuation correction of Ana ~ 1 mag? Why are they so 
similar to those derived using the K98 FIR-luminosity for- 
mula for optically thick starburst galaxies? The distribution 
of Ha attenuations that we derive for the SAPM galaxies is 
similar to that found in other samples of nearby star-forming 
galaxies. To illustrate this. Fig. 5a compares the distribu- 
tion of Aua with those derived by Sullivan et al. (2000) and 
Bell & Kennicutt (2001) using Ha, H/3 and thermal radio 
continuum observations of nearby, UV-selected star-forming 
galaxies [we have converted Ay from figure 1 of Sullivan et 
al. (2000) into Aua = 0.78Av, consistent with the recipe 
these authors used in the first place to compute Av based 
on the Ha/H/3 ratio]. The median Ha attenuation in our 
sample, Aua ~ 1.1, is only slightly larger than those of the 
Sullivan et al. (0.8) and Bell & Kennicutt (0.9) UV-selected 
samples. 

It is important to distinguish here between Ha and V- 



band attenuation. Although most galaxies in our sample are 
optically thick at Ha, they are generally not in the V band. 
This is because the l/-band continuum emission is domi- 
nated by stars that live longer, and hence, that are less ob- 
scured, than those dominating the Ha line emission (Section 
4.1). Fig. 5b shows the distribution of l/-band attenuations 
for our sample. This has a median Av ~ 0.8. The results 
of Fig. 4c then suggest that dust heating by long-lived stars 
roughly compensates the relatively low dust optical depths 
of the galaxies, in such a way that the ratio of FIR luminos- 
ity to SFR is similar to that predicted by the K98 formula 
for optically thick starburst galaxies. 

We now show that differences between the CLOl and the 
K98 Ha estimates of the SFR arise as a result not of one, but 
of several simplifying approximations in the K98 formula. It 
is useful to step back and compare the CLOl SFRs with those 
derived from Equation (|l|), before applying any correction 
for Ha attenuation. This is shown in Fig. 6a. Rather than 
applying a standard attenuation correction of 1 mag, we now 
correct the Ha-derived SFRs of individual galaxies using our 
estimated attenuations shown in Fig. 5 above. As expected, 
this reduces the scatter in the diagram, but an offset of a 
factor of ~ 2 remains between the CLOl and K98 estimates 
(Fig. 6b). A similar discrepancy would presumably be found 
if attenuation were derived, for example, using the Ha/H/3 
ratio. 

An effect that is included in the CLOl model, but not 
in the K98 formula, is the absorption of ionizing photons by 
dust inside H ll regions (Section 4.1 above). We find that, for 
the 149 galaxies in Fig. 6, the mean and rms scatter of the 
fraction of ionizing photons absorbed by dust correspond to 
fhyc ~ 22 ±20 per cent. This is constrained primarily by the 
[S n]/Ha and [O n]/[S ii] ratios of the galaxies, which control 
the dust-to-heavy element ratio and ionization parameter 
of the gas (Section 4.2 and CLOl). For reference, Degioia- 
Eastwood (1992) finds f£^l 35 ± 11 per cent based on 
optical, far- infrared and radio observations of 6 H ll regions 
in the Large Magellanic Cloud. In Fig. 6c, we show the effect 
of correcting the Ha-derived SFRs of individual galaxies in 
Fig. 6b for the fraction of ionizing photons lost to dust. Both 
the scatter and the offset relative to the CLOl estimates of 
the SFR are reduced, but there remains an offset of a factor 
of ~ 1.5. 

Another effect that is included in the CLOl model, but 
ignored by the K98 formula, is the contamination of Ha 
emission by stellar absorption. For the galaxies in Fig. 6, the 
mean and rms scatter of the equivalent width of H-Balmer 
absorption correspond to W^g^= ^ W^f ^ -6.5±1.3 A. This 
compares well with the rough estimate by Kennicutt (1992b) 
of W-i^ ~ —5 A for a sample of 90 nearby spiral galaxies. 
The correction for stellar absorption is therefore quite im- 
portant for galaxies with small observed Ha emission equiv- 
alent widths. This is shown in Fig. 6d, where we correct the 
Ha-derived SFRs of individual galaxies in Fig. 6c for stel- 
lar absorption. This additional correction brings the Ha- 
derived SFRs into agreement with the CLOl model. Some 
minor scatter remains (~ 15 per cent), mainly because the 
CLOl model also accounts for the decrease in the ionization 
rate for increasing stellar metallicity (stars are assumed to 
have the same metallicity as the gas; see Section 4.1 above). 
The K98 prescription, on the other hand, assumes fixed so- 
lar metallicity. We conclude from Fig. 6, therefore, that the 
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neglect of several important physical effects in the derivation 
of the K98 Ha SFR estimator leads to significant underesti- 
mates of the true rate of star formation occurring in typical 
galaxies. 

Fig. 7 presents the estimated SFRs and gas-phase oxy- 
gen abundances of the 705 SAPM galaxies with fluxes and 
equivalent widths of Ha, [O ll], and one or both of [N ll] and 
[Sll]. The SFR in Fig. 7a is plotted per unit bj luminos- 
ity. The median error in SFR/L(fe,j) is about 80 per cent, 
consistent with the results of Fig. 3. Unlike the Ha and 
[On] equivalent widths (Fig. 8), the quantity SFR/L(fej) 
exhibits no dependence on galaxy absolute bj luminosity. 
Many previous studies have noted that galaxies with strong 
Ha or [O ll] equivalent widths dominate the faint end of 
the galaxy luminosity function (e.g., Zucca et al 1997; Love- 
day, Tresse & Maddox 1999; Bromley et al 1998; Blanton & 
Lin 2000; Christlein 2000; Balogh et al. 2001). Our resuhs 
show that this effect must arise in part because the emis- 
sion lines in low-luminosity galaxies are less attenuated by 
dust, rather than because faint emission-line galaxies form 
stars at higher relative rates than their more luminous coun- 
terparts. The gas-phase oxygen abundances of the galaxies 
in our sample increase from 12 -I- log(0/H) « 8.5 ± 0.3 at 
M{bj) ~ -18 to 12 + log(0/H) ^ 8.8 ± 0.2 (i.e. solar abun- 
dance) at M{bj) ~ —21, indicating that more luminous 
galaxies are more metal-rich, although with a large scat- 
ter. This is agreement with other studies of trends in galaxy 
metallicity from observations of H ll regions in nearby star- 
forming systems (e.g. Zaritsky, Kennicutt & Huchra 1994). 



6 SUMMARY AND DISCUSSION 

We have compared a number of different methods for esti- 
mating the star formation rates of galaxies using a variety 
of spectral indicators. Our primary sample is drawn from 
the Stromlo-APM redshift survey and consists of 149 galax- 
ies with fluxes and equivalent widths of Ha, [O ll], and one 
or both of [Nil] and [Sll], as well as IRAS 60 and 100 ^Im 
fluxes. 

We first compared the star formation rates obtained 
using standard estimators based on the Ha, [Oil] and FIR 
luminosities of the galaxies. These estimators were derived 
in inconsistent limits using incomplete models of the repro- 
cessing of radiation from young stars, but they have, never- 
theless, been used by many authors to study star formation 
as a function of galaxy luminosity, morphological type, envi- 
ronment and redshift. Not surprisingly, the three estimators 
give discrepant results. The star formation rates estimated 
using the far-IR fluxes are a factor of ~ 3 higher than those 
estimated using Ha or [On], even after the usual mean at- 
tenuation correction of Aua ~ 1 mag is applied to the data. 

We then derived star formation rates for the same galax- 
ies using estimators constructed from new models developed 
by CLOl, which include a physically consistent treatment of 
the effects of stars, gas and dust. The CLOl star formation 
rates agree much better with those derived from the FIR 
ffuxes than with those derived from the Ha or [O ll] fluxes. 
We showed that the discrepancy between the CLOl and Ha 
estimates of star formation arises primarily because the stan- 
dard Ha estimator does not include the absorption of ion- 
izing photons by dust in H ll regions and the contamination 



of Ha emission by stellar absorption. The strengths of these 
two effects as estimated by the models are in good agree- 
ment with independent observations of normal, nearby star- 
forming galaxies. The agreement with the FIR estimates is 
at least in part a fluke; in this case the two major deviations 
from the assumed opaque-starburst limit have opposite sign 
and approximately cancel for our SAPM sample. 

Our results have important implications for estimates 
of the total integrated star formation density and its evo- 
lution with redshift (Madau et al 1996). Most estimates of 
star formation rate densities at low redshift are based on Ha 
measurements (Gallego et al. 1995; Tresse & Maddox 1998). 
If these studies underestimate the total amount of star for- 
mation in the local Universe by factors of ~ 3, the resulting 
increase of the SFR density to high redshifts, which is gen- 
erally inferred using different estimators, could be consider- 
ably less than commonly believed (see also Cowie, Songaila 
& Barger 1999). We plan to investigate this in more detail 
in a future paper. 

Finally, we have applied the CLOl models to 705 SAPM 
galaxies with fluxes and equivalent widths of Ha, [O ll], and 
one or both of [N ll] and [S ll] to study trends in star forma- 
tion rate and in metallicity as a function of galaxy absolute 
bj magnitude. For this sample, we find that the star for- 
mation rate per unit luminosity is independent of absolute 
magnitude. The gas-phase oxygen abundance does increase 
with galaxy absolute bj luminosity, although the scatter in 
metallicity at fixed magnitude is large. We find that lumi- 
nous galaxies have abundances close to solar on average. 

We note that the star formation rate and metallic- 
ity estimates derived using the CLOl models are based on 
emission-line ratios. In order to obtain reliable estimates of 
these quantities, it is critical that the spectra be accurately 
flux-calibrated. The large wavelength coverage and the ex- 
pected spectro-photometric accuracy of the galaxy spectra 
in the Sloan Digital Sky Survey (York et al 2000) and the 
Two Degree Field survey (Folkes et al 1999) are ideally 
suited to this kind of analysis. These surveys include mea- 
surements of key emission lines not available for the SAPM 
(e.g. [Olll] and H/3), which will greatly strengthen the con- 
straints on the effective ionization parameter and the effec- 
tive dust absorption optical depth, even for galaxies that do 
not have measured far-infrared fluxes. In addition, it will be 
possible to test the CLOl models in detail using the addi- 
tional spectral information. 

In summary, both the models and the data are now 
available that will allow detailed quantitative studies of the 
physical properties of galaxies as a function of luminosity, 
type and environment. This can only lead to an improved 
understanding of the processes that were important in the 
formation of the galaxies we see today. 
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Figure 1. Logarithmic difference between SFRhh, SFR[oii] ^^'^ SFRfir [equations (|lj), 1^ and ^] plotted against absolute rest-frame 
bj magnitude for the 149 SAPM galaxies with fluxes and equivalent widths of Ha, [On], and one or both of [Nil] and [Sll], as well as 



IRAS 60 and 100 /xm fluxes. We have applied a standard attenuation correction of Ajja = 1 mag to SFRhq and SFRi 



the median Alog(SFR) offset for galaxies with M(hj) < 
parentheses). 
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Figure 2. Logarithmic difference between SFRuVi SFRjja and SFRpiR^ [equations (|l|), (|3|) and (^] plotted against absolute rest-frame 
B magnitude for the galaxies with Hq, UV (1600 A) and IRAS 60 and 100 lira flux measurements in the samples of Bell & Kennicutt 
(2001; solid squares) and Meurer et al. (1999; open squares). In each panel, the median Alog(SFR) offset for galaxies with M(bj) < —19 
is indicated along with the associated mean offset and the rms scatter (in parentheses). 
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Figure 3. Logarithmic difference between the CLOl SFR estimates and (a) SFRhcii (b) SFRjon] and (c) SFRpii^ [equations (|l|), 
and (^] plotted against absolute rest-frame hj magnitude for the 149 SAPM galaxies with fluxes and equivalent widths of Ha, [O ll], and 
one or both of [Nil] and [Sll], as well as IRAS 60 and 100 /^m fluxes. The CLOl models are required to fit the observed luminosities and 
equivalent widths of the optical emission lines and the absolute bj magnitude. The filled squares represent galaxies with measurements of 
all four optical lines, the open circles galaxies with Ha, [On] and [Nil] measurements, and the open squares galaxies with Ha, [On] and 
[Sll] measurements. The error bars indicate the uncertainties in the SFRs derived using the CLOl models, as discussed in Section 4.2. 
We have applied a standard attenuation correction of Ana = 1 mag to SFRh^ and SFR[oii]- In each panel, the median Alog(SFR) 
offset for galaxies with M{bj) < — 19 is indicated along with the associated mean offset and the rms scatter (in parentheses). 
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Figure 4. Same as Fig. 3, except tliat the CLOl models are now required to fit the total FIR luminosity (as defined in Section 3), in 
addition to the observed luminosities and equivalent widths of the optical emission lines and the absolute bj magnitude. 
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Figure 5. (a) The distribution of Ha attenuations derived for the 149 SAPM galaxies with fluxes and equivalent widths of Ha, [On], 
and one or both of [Nil] and [Sll], as well as IRAS 60 and 100 /^m fluxes (plain histogram). For comparison, results from Sullivan et 
al (2000) and from Bell & Kennicutt (2001) are shown as shaded and hatched histograms, respectively. Arrows indicate the median 
attenuations in the Sullivan et al. (2000, Aho = 0.8), the Bell & Kennicutt (2001; Ahq = 0.9) and the SAPM (Ahc« = 1-1) samples, (b) 
The distribution of V-band attenuations derived for the same sample of 149 SAPM galaxies as in (a). An arrow indicates the median 
attenuation Ay = 0.8. 
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Figure 6. Breakdown of the difference between the CLOl SFR estimates and SFRhq plotted against absolute rest-frame bj magnitude 
for the 149 SAPM galaxies with fluxes and equivalent widths of Hq, [On], and one or both of [Nil] and [Sll], as well as IRAS 60 and 
100 /tm fluxes. Symbols are as described in Fig. 3. (a) Comparison of the CLOl SFRs with those derived from Equation ), before applying 
any correction for Ha attenuation, (b) After correcting the Ha-derived SFRs of individual galaxies using the estimated attenuations 
shown in Fig. 5. (c) After further correction for the fraction of ionizing photons absorbed by dust in Hll regions, (d) After further 
correction for stellar absorption. In each panel, the median Alog(SFR) offset for galaxies with M{bj) < —19 is indicated along with the 
associated mean offset and the rms scatter (in parentheses). 
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Figure 7. Estimated SFRs and gas-pliase oxygen abundances of the 705 SAPM galaxies with fluxes and equivalent widths of Ha, 
[On], and one or both of [Nil] and [Sll] plotted against absolute rest- frame bj magnitude. Symbols are as described in Fig. 3. (a) Star 
formation rate per unit absolute bj luminosity, (b) Gas-phase oxygen abundance [12 + log(0/H)Q ^ 8.8]. For clarity, only 1 in 7 error 
bars are shown to indicate the uncertainties in the SFRs and gas-phase oxygen abundances derived using the CLOl models (as discussed 
in Section 4.2). 
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Figure 8. Ha and [O ll] rest-frame emission equivalent widths of the 705 SAPM galaxies with fluxes and equivalent widths of Ha, [O ll] , 
and one or both of [N ll] and [S ll] plotted against absolute rest-frame bj magnitude. For clarity, only 1 in 7 observational error bars arc 
shown. 
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